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Introduction: The small obliquity of Mercury causes to- 
pographic depressions located near its poles to cast persis- 
tent shadows. Many [1, 9, 15] have shown these perma- 
nently shadowed regions (PSRs) may trap water ice for ge- 
ologic time periods inside cold-traps. More recently, direct 
evidence for the presence of water ice deposits inside craters 
was remotely sensed in RADAR [5] and visible imagery [3]. 
Albedo measurements (reflectence at 1064 nm) obtained by 
the MErcury Space ENviroment GEochemistry and Ranging 
Laser Altimeter (MLA) found unusually bright and dark ar- 
eas next to Mercury’s north pole [7]. Using a thermal illumi- 
nation model, Paige et al. [8] found the bright deposits are 
correlated with surface cold-traps, and the dark deposits are 
correlated with subsurface cold-traps. They suggested these 
anomalous deposits were brought to the surface by comets 
and were processed by the magnetospheric radiation flux, re- 
moving hydrogen and mixing C-N-O-S atoms to form a vari- 
ety of molecules which will darken with time. Here we use a 
thermal illumination model to find the link between the cold- 
trap area fraction of a rough surface and its albedo. Using this 
link and the measurements obtained by MESSENGER we de- 
rive a surface and a subsurface ice distribution map on Mer- 
cury’s north pole below the MESSENGER spatial resolution, 
~ 500 m. We find a large fraction of the polar ice on Mer- 
cury resides inside micro cold-traps (of scales 10 — 100 m) 
distributed along the inter-crater terrain. 


The Bimodal Albedo Distribution: When we examine 
the albedo distribution of Mercury’s north pole we find it is 
bimodal with an extended tail (Figure 1). The dark regions 
have an albedo < 6% while the bright regions have an albedo 
> 25%. We use this distribution to model the surface re- 
flectance, which will be explained next. To describe the mean 
albedo of the surface we use the mean of the higher mode, 
~ 16%, and to describe the anomalous dark deposits we use 
the mean of the lower mode ~ 4%. We consider slopes as 
bright if their albedo is higher than 25%. 


Modeling Rough Surfaces: In the absence of atmo- 
spheric erosion, topography on airless bodies can be divided 
into craters and inter-crater terrain. While craters have a well 
defined, closed shape (in map view), inter-crater topography 
is more random in nature. A common method in which this 
type of topography can be modeled is by using random Gaus- 
sian surfaces [2, 6, 14]. Such surfaces are created by mod- 
ifying the power spectrum of a random field (2D) of nor- 
mally distributed slopes according to measured parameters 
such as the Hurst exponent [13] and RMS slope o,. Since 
the Hurst exponent and RMS slope are unknown for Mercury 
at our scales of interest, we assume values measured on the 
Moon [10]; surfaces in scales 10 — 100 m have an RMS slope 
a, = 10°—15°. An example for such a surface with 0, = 15° 
can be seen in Figure 2. 
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Figure 1: The albedo (reflectance at 1064 nm) distribution 
of Mercury’s north pole between latitudes 75° and 84°. The 
histogram shows the probability density function (PDF) with 
a bin width of 0.01. 
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Figure 2: An example for the modeled maximum tempera- 
tures on a modeled rough surface with 0, = 15°. 


Modeling Cold-Traps: We choose to define a cold- 
trap as an area in which the mass loss rate of water ice 
< 1cm Ga. This mass loss rate (and thus the stability of 
cold-traps) is extremely sensitive to changes in temperature. 
For example, for the mass loss rate we chose, a surface cold- 
trap must not exceed a maximum temperature of 110K [15]. 
Subsurface cold-traps are more stable than surface cold-traps 
due to the diffusion barrier created by the overlaying regolith, 
and are modeled by considering the dependence of the mass 
loss rate on depth and temperature [12]. To find the maxi- 
mum temperature distribution on and below a rough surface, 
we employ a thermal illumination model accounting for in- 
solation, scattering, thermal emission and conduction into the 
subsurface in 1D [11]. The thermal parameters we use (con- 
ductivity, heat capacity) vary with depth and temperature, but 
the typical thermal inertia and skin-depth are ~ 25SI and 
~ 10 cm. Using the derived maximum temperature distri- 
bution, we find the cold-trap area fraction on and below the 
surface. Next, we link the cold-traps area fraction and the 
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Figure 3: Comparison between the MESSENGER measured 
mean albedo (over a latitude ring) and the results of the model 
(LHS y-axis) assuming two roughness values, 0, = 10° and 
go; = 15°. The RHS y-axis shows the corresponding subsur- 
face cold-trap area fraction. 


albedo on rough surfaces. 

Linking Cold-Traps and Albedo: To model the surface 
reflectivity we assign surface cold-traps with an albedo of 
25%, subsurface cold-traps with an albedo of 4% and non 
cold-trapping regolith with an albedo of 16%. The over- 
all surface albedo will simply be the linear combination of 
the albedo of the slopes composing it [4]. In Figure 3 we 
show the modeled albedo compared to the latitudinal mean 
albedo measured by MESSENGER. As we do not account for 
craters, our modeled albedo slightly deviates from the mean 
measured albedo. However, our model predicts the 2% de- 
crease in albedo between latitude 75° — 84° witnessed by 
MESSENGER. On the RHS y-axis (red) we show the cor- 
responding area fraction of subsurface cold-traps. We find 
that in polar latitudes, rough random Gaussian surfaces have 
a subsurface cold-trap area fraction of up to ~ 20%. 

Micro cold-traps on Mercury’s north pole: Finally, we 
employ the correlation we found above and the results of a 
thermal illumination model at larger scales to create maps 
of surface and subsurface ice below MESSENGER’s spatial 
resolution. Figure 4 is an example for one of the maps we 
created, showing the distribution of subsurface ice on Mer- 
cury’s poles - assuming the roughness in scales of ~ 10 m is 
oa; = 10°. As previously shown [8], we find ice is stable on 
the permanently shadowed, pole-facing slopes of craters due 
to their protective closed shape. We find subsurface ice may 
persist even in areas outside craters, where the maximum tem- 
perature distribution is controlled by the small scale rough- 
ness and PSR distribution. For a given latitude, the number of 
permanently shadowed slopes increases when the roughness 
increases and slopes with surface cold-traps replace slopes 
with subsurface cold-traps. However, the number of subsur- 
face cold-traps increases as well, and the overall albedo of the 
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Figure 4: One of the maps produces using our derived link be- 
tween albedo and ice stability; the distribution of subsurface 
ice on Mercury’s pole in area smaller than the MESSENGER 
scale, assuming the roughness in scales 10m is 0, = 10°. 


surface slightly decreases. Assuming Mercury’s polar inter- 
crater terrain is a Gaussian random surface with RMS slope 
os = 10°, we find it holds ~ 3% subsurface ice (below dark, 
weathered regolith) and ~ 1% surface ice. We conclude a 
significant amount of the ice on Mercury’s north pole is con- 
tained inside micro cold-traps created by small-scale rough 
features. While most of the polar ice is trapped on the per- 
manently shadowed pole-facing slopes of craters, ice within 
small scale (micro) cold-traps is potentially more accessible 
and therefore easier to investigate. 
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